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(57) Abstract 

A transmit modulator which uses a quadrature modulator (102) has a predistortion block (100) which produces a predistorted 
output signal. The predistortion block (100) pre-compensates for errors introduced by the quadrature modulator (102) based on a set of 
predistortion i parameters. The quadrature modulator (102) receives the output signals from the predistortion block (100). The quadrature 
modulator (102) up~converts the I and Q channel signals and combines them. In the process, the quadrature modulator (102) introduces 
erT £!" S ; • , caIlbratlon mode > a sinusoidal wave at frequency f ca , is applied to the input of the predistortion block (100) (while the modulator 
output is downconverted with an AM-detector). Spurious energy produced by the quadrature modulator errors but reduced by the effect 
of the predistortion block (100) is generated at f^i and 2*f ca i. The quadratic polynomial minimization calculator (166) determines a 
subsequent value of the predistortion parameters based upon a quadratic relationship between energies present in the digital representation 
of the spectrum of the output of the quadrature modulator (102) at frequencies f ca i and 2*f ca i and the values of the previous predistortion 
parameters. 
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METHOD AND APPARATUS FOR THE DETERMINATON OF THE PREDISTORTION PARAMETERS IN A QUADRA- 
TURE MODULATOR, BY MEASURING THE OUTPUT POWER 

Background of the Invention 

I. Field of the Invention 

This invention relates generally to wireless transmitters. More specifically, the invention relates to 
quadrature modulation of a transmission signal in a wireless transmitter. 

II. Description of the Related Art 

Modern digital communications systems often use quadrature techniques to impress data on a modulated 
signal. The use of quadrature techniques allows twice as much data to be transferred within the same bandwidth 
compared to single-phase modulation techniques. Quadrature techniques can be applied in many types of systems 
using a variety of modulation and access techniques. For example, quadrature techniques can be applied in code 
division multiple access (CDMA), time division multiple access (TDMA) and frequency division multiple access (FDMA) 
systems as well as others. Using quadrature techniques, at the transmitter, a quadrature modulator impresses a 
portion of the signal energy on an in-phase (I) channel and the remainder on a quadrature (Q) channel which is 90 
degrees out of phase in comparison with the I channel. At the receiver, the energy from each channel can be 
separately recovered by a quadrature demodulator. 

In many modern wireless communication systems, the cost of the remote communication unit is a 
substantial barrier to extensive deployment. For example, in a satellite based, wireless local loop telephone system, 
the cost of the handset to the end consumer can be a major factor in determining the penetration into the market 
which the wireless service will enjoy. For this reason, considerable effort to techniques has been expended in recent 
years revolutionize remote unit architectures so that they cost less. 

In a typical wireless transmitter, a digital representation of a baseband signal is generated by digital 
circuitry. An analog to digital converter (A/D) is used to convert the digital representation to a baseband analog signal. 
The baseband signal is upconverted from baseband to a fixed intermediate frequency (IF) where a variety of signal 
processing functions are applied. For example, the gain of the signal can be set to accommodate current system 
conditions. The signal can be subjected to rigorous filtering in order to reduce transmission spurs. The conditioned IF 
signal is converted to a channel sensitive transmission radio frequency (RF). The RF signal is transmitted over the 
wireless link. 

The main issues in designing a quadrature modulator are maintenance of quadrature phase between the I and 
Q channel, minimization of differential gain errors and differential direct current (D.C.) offsets between the I and Q 
channels as well as minimization of carrier leakage. Each of these factors can cause spurious output power to be 
generated which can interfere with the transmitted signal as well as signals in adjacent bands. In addition, each of 
these factors can introduce distortion into the transmitted signal which can result in an increased bit-error rate (BER) 
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at the receiver and an undesired spreading of the signal spectrum. In practical circuit designs, these factors vary with 
temperature, device biasing, component aging and frequency, often making readjustment during operation necessary. 
In addition to baseband signal processing, IF signal processing can be used to accommodate for some of these errors. 

In order to reduce the cost of the remote unit, the luxury of using IF techniques is no longer economically 
practical in many situations. Alternatively, modern communication transmitters have begun to use direct conversion 
techniques which provide for the direct conversion of a baseband signal to a radio frequency signal. Direct conversion 
systems are characterized by the use of a single up-conversion stage using a high-frequency quadrature modulator. 

To construct a high frequency, direct conversion circuit, it is necessary to mount many extremely small 
devices on a single semiconductor chip substrate. These devices are closely packed together and, further, because of 
their size, are extremely sensitive to stray currents. Because of the dense packing and the associated network of 
interlaced wires, the circuits are susceptible to crosstalk between the various components. This crosstalk is primarily 
due to capacitive coupling between adjacent wires, but may also result from inductive coupling and transmission line 
effects. In addition, D.C. offsets and gains realized even on adjacent signal paths on a single substrate can vary 
significantly. Although through careful circuit design and layout these factors can be controlled, the unpredictability 
of extremely complex systems results in random variations in these factors from path to path, from part to part and 
from assembled system to assembled system. 

Despite these design challenges, direct conversion architectures provide several desirable features. For 
example, direct conversion techniques tend to require less circuitry leading to a higher efficiency and reduced D.C. 
power requirements. Direct conversion transmitters are typically less costly to manufacture and smaller in size than 
conventional systems. Figure 1A is a block diagram illustrating a basic direct conversion transmitter. As shown, a 
transmit modem (TM) 10 feeds a quadrature modulator (QM) 14 which in turn drives a power amplifier (PA) 16. The 
transmit modem 10 generates the complex baseband information signals. The quadrature modulator 14 provides 
direct conversion of the combined complex baseband information signals to the RF transmit frequency. The power 
amplifier 16 amplifies the RF signal for transmission over the wireless link. 

Because a direct conversion system is very sensitive to quadrature modulator gain imbalance, phase 
imbalance and D.C. offsets, careful control over quadrature modulator errors is required. If quadrature modulator 
errors were more or less static, then, a simple calibration of the phase, gain, and offset errors would suffice. 
However, as noted above, the quadrature modulator errors can be expected to change with temperature, channel 
frequency, device biasing and component aging, and, thus, some means of "tracking" the quadrature modulator errors 
is needed. 

Several error correction techniques have been studied in the literature such as Cartesian feedback, 
feedforward and predistortion. Using predistortion, the baseband signal is passed through a nonlinear system having 
an inverse characteristic compared to those of the quadrature modulator causing the overall system to be linear. 

Figure 1B is a block diagram illustrating a basic direct conversion transmitter which incorporates 
predistortion. In comparison with Figure 1A, a predistorter 12 has been inserted between the transmit modem 10 and 
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the quadrature modulator 14. Ideally, the predistorter 12 exhibits the inverse transfer function of the undesired 
characteristics of the quadrature modulator 14. The transmit modem 10 generates the complex baseband information 
signal which is then pre-corrected in the predistorter 12 to compensate for anticipated errors which will be introduced 
by the quadrature modulator 14. 

5 Several approaches have been used to determine the value of the transfer function of the predistorter. In 

JX Cavers, "A Linearizing Predistorter with Fast Adaptation," IEEE Trans. Vehicular Technology, vol. 39, no. 4, pp. 
374-382, Nov. 1990, a look-up table having input power as the index and complex gain as the table entries is used. 
The table look-up method can be used to fit any gain profile, and performance can be adjusted by increasing or 
decreasing the table size. The second approach, as detailed in S.P. Stapleton, G.S. Kandola, and J.K. Cavers, 
1 0 "Stimulation and Analysis of an Adaptive Predistorter Utilizing a Complex Spectral Convolution," IEEE Trans. Vehicular 

Technology, vol. 41, no. 4, pp. 387-394, Nov. 1992., uses two polynomial equations to fit the desired complex gain 
curves. This method is simpler to implement, but arbitrary gain profiles may be difficult to fit with polynomial 
functions. The polynomial method, because it has fewer variables, is also easier to initialize and update than the table 
look-up method. 

15 A bin a r y search is proposed by R. Datta and S.N. Crozier, "Direct Modulation at L-Band using a Quadrature 

Modulator with Feedback," International Mobile Satellite Conference, Ottawa (IMSC '95) pp. 383 - 388. The binary 
search is used to find the predistortion parameters by the measurement of energies at two or more frequencies during 
the calibration process. The binary search method is an iterative process which continues until the process converges 
on an acceptable value. The iterative nature of the binary search process can consume a significant amount of system 
20 resources and require significant amount of time. In addition, the binary search algorithm converges to the correct 

solution only if the correct solution is within the initial range of uncertainty. 

Each of the prior art methods requires a substantial amount of processing power and can be exceedingly 
difficult to implement. In addition, they may each require a substantial amount of time to converge upon the optimal 
values and may not be robust under certain operating conditions. 
25 Therefore, there has been a need in the industry, to develop a robust and reliable determination of 

predistortion parameters which does not consume an excessive amount of processing power. 

Summary of the Invention 

A transmit modulator which uses a quadrature modulator has a predistortion block which produces a 
predistorted output signal. The predistortion block pre-compensates for errors introduced by the quadrature 
30 modulator. The predistortion block comprises a D.C. offset compensation summer which determines a D.C. offset 

adjustment to an I channel signal path based upon a predistortion parameter, bj. The predistortion block also has a 
D.C. offset compensation summer which determines a D.C. offset adjustment to a Q channel signal path based upon a 
predistortion parameter, b q . The predistortion block also has a crosstalk amplifier which determines a crosstalk gain 
between the I channel and the Q channel signal path based upon a predistortion parameter, K. The predistortion block 
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The quadrature modulator receives the output signals from the predistortion block. The quadrature 
modulator up-converts the I and Q channel signals and combines them. In the process, the quadrature modulator 
introduces errors. These errors create spurious energy. 

In calibration mode, a sinusoidal wave at frequency f ca , is applied to the input of the predistortion block and a 
transformer is coupled to the output of the quadrature modulator. The transformer produces a digital representation 
of a spectrum of the output of the quadrature modulator converted to baseband. Spurious energy produced by the 
quadrature modulator errors but reduced by the effect of the predistortion block is generated at f ca , and 2*f cal . 

A quadratic polynomial minimization calculator receives the output of the transformer. The quadratic 
polynomial minimization calculator determines a subsequent value of the predistortion parameters, b lt b q , K and G 
based upon a quadratic relationship between energies present in the digital representation of the spectrum of the 
output of the quadrature modulator at frequencies f ca , and 2*f cal and the values of the predistortion parameter b i( b q , K 
andG. 

Brief Description of the Drawings 

The features, objectives, and advantages of the invention will become more apparent from the detailed 
description set forth below when taken in conjunction with the drawings wherein like parts are identified with like 
reference numerals throughout and wherein: 

Figure 1A is a block diagram illustrating a basic direct conversion transmitter. Figure 1B is a block 

diagram illustrating a basic direct conversion transmitter which incorporates predistortion. 

Figure 2 is a block diagram showing a predistortion quadrature modulator incorporating the invention. 

Figure 3 is a flowchart showing the iterative process by which a set of predistortion parameters is selected. 

Detailed Description of the Invention 

It is known from basic quadrature modulator theory that output spurs are created by the error introduced by 
a quadrature modulator. It is also known that the energy at the spur frequencies is essentially defined as a quadratic 
function of those errors. For this reason, the energy at the spur frequencies is also defined by a quadratic function of 
any predistortion parameters which are applied before quadrature modulation. The prior art does not exploit this 
knowledge in the predistortion parameter determination process. The invention uses this knowledge to develop an 
efficient, iterative process which is easy to implement and which converges reliably under many conditions. The 
invention fits a quadratic polynomial to the measured energy at several spur frequencies and, then, finds the abscissa 
value at which this polynomial assumes an extremum. From theory, we know that this extremum must be a minimum. 
This process can be executed iteratively until sufficient convergence is obtained. 

Figure 2 is a block diagram showing one embodiment of a predistortion quadrature modulator incorporating 
the invention. The predistortion quadrature modulator shown in Figure 2 performs direct modulation of a signal at the 



WO 00/31941 



PCT/US99/27256 



frequency of transmission. The predistortion quadrature modulator of Figure 2 operates in one of two modes. The 
first mode is a transmission mode in which the signals to be transmitted are modulated according to normal operation 
and transmitted over a wireless link. In transmission mode, the circuit uses the previously determined predistortion 
parameters. The second mode is a calibration mode in which no signal is being transmitted. During calibration mode, 
the predistortion parameters are determined. In a TDMA system, calibration mode can be active between the 
transmission slots. In a FDMA, CDMA or TDMA system, calibration mode can be active between transmission bursts. 

The predistortion quadrature modulator shown in Figure 2 is comprised of four blocks. A digital predistortion 
block 100 pre-compensates the baseband transmission or calibration mode signals and converts them to analog form. 
A quadrature modulation block 102 represents the quadrature modulator including the imperfections which it 
introduces. In calibration mode, an amplitude modulation (AM) detection block 104 uses an AM detection process to 
generate a baseband analog signal. In transmission mode, the AM detection block 104 simply passes the transmission 
signal to the next stage. In calibration mode, a parameter determination block 106 collects a series of data samples 
and executes a series of calculations in order to determine a set of predistortion parameters for use in subsequent 
transmission mode operation. 

In transmission mode, the input to the digital predistortion block 100 is the I channel and Q channel signals 
in baseband digitized form. In calibration mode, the input to the digital predistortion block 100 is a digitized sinusoidal 
wave at a frequency f Ml . The digital predistortion block 100 contains four adjustable elements which are used to pre- 
compensate for errors introduced by circuit imperfections as representatively shown in the quadrature modulation 
block 102 in Figure 2. During transmission mode, the adjustable elements are set according to the predistortion 
parameters most recently determined. In calibration mode, the settings of the four adjustable elements are varied over 
a series of values in order to determine a new set of predistortion parameters for use during a subsequent 
transmission mode. 

The Q channel input is coupled to the input of a crosstalk amplifier 1 1 2. The I channel input is coupled to a 
summer 110 which sums the I channel input with the output of the crosstalk amplifier 112. The gain of the crosstalk 
amplifier 1 1 2 is determined by the predistortion parameter K. The cross talk amplifier 1 1 2 couples a small portion of 
the Q channel signal into the I channel signal path to pre-compensate for cross talk imperfections. Crosstalk may also 
be modeled as a phase shift of I channel signal path relative to the Q channel signal path. The crosstalk amplifier 112 
can also be implemented by coupling a portion of the I channel signal into the Q channel path in an alternative 
embodiment. 

The output of the summer 1 10 is coupled to a unity gain amplifier 1 14. The Q channel input is coupled to a 
gain predistortion amplifier 116. The gain of the gain predistortion amplifier 116 is determined by the predistortion 
parameter G,. Because it is only necessary to balance the relative gain between the I channel and the Q channel 
paths, it is only necessary to compensate the gain of one of the signal paths. In this case, the Q channel is arbitrarily 
selected to comprise the gain predistortion amplifier 116. 
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The output of the unity gain amplifier 114 and the gain predistortion amplifier 116 are coupled to a D.C. 
offset predistortion summer 118 and a D.C. offset predistortion summer 120, respectively. The D.C. offset - 
predistortion summer 118 sums a D.C. offset signal determined by the predistortion parameter b. The D.C. offset 
predistortion summer 120 sums a D.C. offset signal determined by the predistortion parameter b q . The output of the 
5 D.C. offset gain predistortion summers 118 and 120 are coupled to digital-to-analog (D/A) converters 122 and 124, 

respectively. The D/A converters 122 and 124 produce low frequency or baseband analog signals which are coupled 
to the quadrature modulation block 102. 

The quadrature modulation block 102 is comprised of several actual circuit elements as well as a set of 
representative elements which represent the errors introduced within a quadrature modulator. These representative 
10 elements illustrate the introduction of undesirable offsets, ingresses and imbalances due to circuit imperfections. 

These representative elements are not actual circuit elements which would be found in an actual implementation, but 
are, instead, representations of the undesirable effects added to the figure in order to facilitate an understanding of 
the predistortion mechanism. In the quadrature modulation block 102, a representative summer 130, a representative 
summer 132, a representative amplifier 134, a representative amplifier 136 and a representative phase shifter 140 
1 5 each illustrate an error introduced by the quadrature modulator 1 02 for which predistortion is included. None of these 

elements is an actual circuit element which would be included in an actual implementation. In an actual circuit 
implementation, the representative elements are replaced with short circuits. 

The imperfections modeled by the representative elements as well as the other circuit elements in the 
quadrature modulation block 102 operate the same way in transmit mode as they do in calibration mode. The 
20 imperfections modeled by the representative elements may vary as a function of time due to temperature changes and 

other environmental changes as well as aging. For this reason, the calibration process continues to be repeatedly 
executed during operation of the circuit. 

The outputs of the D/A converters 122 and 124 are coupled to the representative summers 130 and 132, 
respectively. The undesirable D.C. offset for the I channel a, is summed with the I channel signal by the representative 
25 summer 130. The undesirable D.C. offset for the 0 channel a, is summed with the 0 channel signal by the 

representative summer 132. The outputs of the representative summers 130 and 132 are input into the 
representative gain amplifiers 134 and 136, respectively. The gain of the representative amplifier 134 is g, The gain 
of the representative amplifier 134 is g,. In general, these gains are not equal to one another, thus, introducing a gain 
imbalance error into the modulator. 
30 The outputs of the representative amplifiers 134 and 136 are coupled in the inputs of mixers 142 and 144, 

respectively. The local oscillator inputs of the mixers 142 and 144 are driven by a local oscillator generation unit 138. 
The actual frequency that the local oscillator produces depends on how the quadrature modulator is used in the actual 
transmitter circuitry. In one embodiment, the frequency of the local oscillator is an intermediate frequency (IF). In 
another embodiment, the frequency of the local oscillator is the actual RF center frequency of transmission. In yet 
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another embodiment, the frequency of the local oscillator is programmable and may be assigned a different frequency 
in calibration mode than in transmit mode. 

The representative phase shifter 140 illustrates the phase imbalance between the I and Q channel and is 
also representative of the effects of crosstalk in the system. The outputs of the mixers 142 and 144 are summed 
together in a summer 148. The output of the quadrature modulation block 102 is a signal at RF or IF which is the sum 
of the I and Q channel signals. 

In transmit mode, a switch 154 within the AM detection block 104 is set to the A position and the combined 
I and Q channel signal is passed to the next transmission stage (not shown). For example, in one embodiment, in 
transmit mode, the AM detection block 104 routes the analog signal produced by the quadrature modulation block 102 
to an upconversion stage or to a power amplifier. 

In calibration mode, the switch 154 is set to the B position so that the output of the quadrature modulation 
block 102 is coupled to an AM detector 156. In one embodiment, the AM detector 156 squares the input signal which 
essentially moves the spectrum of the signal to baseband. The AM detector 156 may be replaced with other well- 
known circuit elements which translate an RF signal to a lower frequency signal. 

In calibration mode, any energy component which appears at the input of the AM detection block 104 at the 
frequency of the local oscillator minus f cal is present at 2*f eal at the input to the parameter determination block 106. 
Any local oscillator leakage component caused by D.C. imbalance which appears at the input of the AM detection 
block 104 at the frequency of the local oscillator is now present at f caI at the input to the parameter determination 
block 106. Any energy component caused by intermodulation nonlinearities which appears at the input of the AM 
detection block 104 at the third harmonic frequency is present at f taf at the input to the parameter determination block 
106. 

Within the parameter determination block 106, a low-pass filter 160 filters the incoming signal. An analog- 
to-digital <A/D) converter 162 converts the analog signal to a digital signal. A Fast Fourier Transform (FFT) 164 
produces a digital spectral representation of the signal. In an alternative embodiment, the FFT 164 can be replaced 
with a Discrete Fourier Transform (DFT) or any other mechanism which produces a digital representation of the 
spectrum of the signal. 

The quadratic polynomial minimization (QPM) calculator 166 executes the QPM algorithm. The QPM 
calculator 166 determines updated values for the predistortion parameters {b„, G v b f and K) used in block 100 during 
transmission mode based upon the data collected during calibration mode. One exemplary embodiment of the 
operation of the QPM calculator 166 is described in detail below with relation to Figure 3. The OPM calculator 166 
outputs the predistortion parameters to the corresponding circuit element in the digital predistortion block 100, thus, 
completing the calibration feedback loop. 

When an installed unit is operating in the field, routine operation is interrupted to enter a calibration mode. 
As noted above, in a TDMA system, calibration mode can be entered between transmission slots. In a TDMA, FDMA 
or CDMA or other type of system, calibration mode can be entered during a natural period of idleness which occurs 
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between the transmission bursts. If the length of the period of idleness does not allow sufficient time for the 
collection of a sufficient amount of data for execution of the calibration process, the collection of data can be 
executed over a series of disjoint periods of calibration. Alternatively, transmission of data can be interrupted or 
suspended in order to accommodate the calibration process. Subsequent to data collection, the actual determination 

5 process may continue as the circuit re enters transmission mode. 

Typically, the QPM calculator 166 is implemented in software running on a digital signal processor (DSP) or 
is embodied in hardware such as in a field programmable gate array (FPGA) or in an application specific integrated 
circuit (ASIC). In one embodiment, the QPM calculator 166 comprises a series of processes that when executed 
perform the steps shown in Figure 3. No matter what media and mode is selected to implement the invention, in order 

1 0 to reduce the cost, the operation of the QPM calculator should be simple and short. The operation detailed in Figure 3 

satisfies these requirements. 

The operation of the QPM calculator 166 is an iterative process although in actual practice the invention 
often converges with one iteration . Each iteration of the process consists of two distinctive substages. Each 
substage requires a set of three measurements of energies, each measurement corresponding to a set of predistortion 
1 5 parameters. Each measurement is implemented by feeding a sinusoidal calibration signal with a frequency of f ca) into 

the digital predistortion block 100 and collecting a set of signals output by the AM detection block 104 over a 
predetermined period of time. In a typical embodiment, the period of time for sample collection is an integer multiple 
of the period of f cal where the integer multiple is on the order of 4 to 10. In one embodiment, the sample sets are 
stored in a buffer for subsequent processing in non-real time. In another embodiment, if there are sufficient central 
20 processing unit resources available, the samples are processed in real time. 

In general, the energy at f^ is affected by changes in b ( and b, while the energy at f — is not significantly 
affected by changes in G„ and K. Likewise, in general, the energy at f 2cal is affected by changes in G, and K while the 
energy at f 2ca , is not significantly affected by changes in b, and b„. Thus, in general, in each set of measurements, one 
of the parameters varies the energy at f cat and one of the parameters varies the energy at f 2cal . Sets of three 
25 calibration measurements are performed. Within each set, one parameter affecting the energy at f cal and one 

parameter affecting the energy at f 2caI are varied. For illustration purposes, the example that follows groups b, with G, 
and K with b. The resulting three sets of measurements for the first and second substage are shown in Table I and 
Table II, respectively. 
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The first time the iterative process is executed, the value of G,(n), b,(n), K(n) and b ( (n) are set to an arbitrary 
value within a suitable range of uncertainty for these parameters. In one embodiment, the value of the predistortion 
parameter G q is initially set to unity and the value of the other parameters are initially set to 0. After an initial 
determination has been made, each time the iterative process is re-executed, the previous values of the predistortion 
parameters are used as a starting point. In a typical system, the changes in the imperfections of the quadrature 
modulator occur at such a slow rate that the previous values of the predistortion parameters are typically a suitable 
starting point for the determination of a new set of predistortion parameters. For each iteration of the first substage, 
step values Ab q and AG, is selected corresponding to b, and G„, respectively. The value of these parameters may 
decrease as the number of iterations increases. A typical value for Ab, is about 20% of the amplitude of the 
calibration signal and a typical value for AG, is about 20% of the nominal gain of the corresponding gain compensation 
amplifier. Likewise, for the second substage, step values Abj and AK must be chosen corresponding to bj and K, 
respectfully. A typical value for AK sets about a 10% crosstalk factor between the two channels (e.g. about 10% of 
the Q channel signal is cross fed into the I channel) and a typical value of Ab, is also about 20% of the amplitude of 
the calibration signal. 

Referring to Figure 3, a description is presented of the iterative process by which a set of predistortion 
parameters is determined based upon the essentially quadratic relationship between the predistortion parameter and 
the spurious energies generated by the imperfect quadrature modulator. Figure 3 is exemplary and other variations 
which related the predistortion parameter to the spurious energies based upon a set of quadratic equations can be 
developed in accordance with these teachings. 

The flowchart of Figure 3 represents operation of the quadrature modulator in calibration mode, in general, 
although some of the steps which do not involve the collection of data may be also executed in transmission mode. 
The flow begins in start block 200. In block 202, a sinusoidal wave at frequency f C3l is applied in digital form to the 
digital predistortion block 100 and a set of samples are collected with the predistortion parameters set at an initial 
value. For example, in one embodiment, the D.C. offset predistortion parameters, b, and b,, are set to zero, the gain 
imbalance predistortion parameter, G, is set to unity, and the crosstalk predistortion parameter, K, is set to zero such 
that no predistortion is provided. The energy output by the FFT 164 at f tal and 2 # f c> , is determined and the values are 
stored for later reference. 

In block 204, samples corresponding to three sets of predistortion parameters shown in Table I are 
collected. The first set of samples is collected with the predistortion parameter G, set equal to the previous value 
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offset on the low side by the selected step value ~G q . For notation purposes, we shall designate this value as x T 2 ' cal . 
The predistortion parameter b q is set equal to the previous value offset on the low side by the selected step value ~b r 
The value of the predistortion parameter K and the predistortion parameter bj are each set equal to their respective 
previous value. For notation purposes, we shall designate this value as x/ cai . 

For the second set of samples, all four predistortion parameters are set equal to their previous value. In one 
embodiment, during the first iteration, this data is the same as that which is collected in block 202 and that data can 
simply be re-used rather than re-derived from samples. For notation purposes, we shall refer to the previous value of 
the predistortion parameter G q (n) as x 2 2fcal and the previous value of the predistortion parameter b q (n) as x 2 fcal . 

The third set of samples is collected with the predistortion parameter G q set equal to the previous value 
offset on the high side by the selected step value ~G q . For notation purposes, we shall refer to this value as x 3 2,cal . 
The predistortion parameter b q is offset on the high side by the selected step value ~b q . For notation purposes, we 
shall refer to this value as x™. The first three sets of samples are used in the first substage of the iterative process. 

In block 206, the energy at f ca , is determined for the first three sets of samples. For notation purposes, we 
shall refer to the energy values as y/ cal , y 2 fcal , and y 3 fcal for the sample sets 1, 2 and 3, respectively. In block 208, 
three quadratic equations are developed as shown in Equation 1 below,, based upon the theory stated above that the 
energy at the spur frequencies is essentially a quadratic function of predistortion parameters. 

y/ ca1 = aix?*) 2 + bx* 31 + c for i - 1, 2 and 3 Equation 1 

Also in block 208, the value.of a and b are determined. For example, the formulas in Equations 2 and 3, below, may 
be used to solve for the value of a and b based upon the three quadratic equations given by Equation 1. 

a «= [y/ ca, (x 2 ,cal -x 3 ,c *) • Yj'^x/'V) + y 3 ,ca, (x/ ca, -x 2 fca1 )] / Det Equation 2 

b- [(x 1 fca, ) 2 (y 2 ,ca, -y 3 ,ca, ) - (x 2 ,ca '} 2 (y/ c V ai ) + <x 3 fcal ) 2 (y/ ca, -y 2 fcal )] / Det Equation 3 

where 

Det - (x/ cal ) 2 (x 2 ,ca, -x 3 fcal ) - (x 2 ,ca! ) 2 (x 1 ,ca, -x 3 ,caI ) + (x 3 ,cal ) 2 (x/ ca, .x 2 fcal ) 

In block 210, the value of the predistortion parameter b q for the next iteration is determined as shown in Equation 4, 
where the minimum value of the polynomial y fcal occurs. 

b q (n+l) - -b/(2a) Equation 4 

Note that the determinant (Det) cancels out of Equation 4. Therefore, the value of the determinant does not have to 
be calculated in an actual implementation. 

One principle upon which the invention operates is that a functional relationship exists between the spurious 
energy at f ca , and the value of the predistortion parameter b q . Stated another way, a value of the predistortion 
parameter b q exists for which the associated energy at f C3) is a minimum. The same is true for predistortion parameter 
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bj. In a similar manner, a value of the predistortion parameters G q , exists for which the associated energy at f 2caI is a 
minimum. The same is true of the predistortion parameter K. 

Still using the first three sets of samples, in block 212, the energy at 2 # f caf is determined. For notation 
purposes, we shall refer to the energy values as y, 2 '"', y 2 2,cal , and y 3 2fcal for the sample sets 1, 2 and 3, respectively. In 
block 214, three quadratic equations are developed as shown in Equation 5. 

y iw . a(x 2fca!,2 + bx 2iw + c for i = 1 , 2 and 3 Equation 5 

In block 216, the value of a and b are determined. For example, the formulas given in Equations 6 and 7 may be used 
to solve for the value of a and b based upon the three quadratic equations given by Equation 5. 

a = [ yi 2fca, (x 2 2fcal .x 3 2,c 1 • y 2 2,ca! (x 1 2^ca, -x 3 2,£a, ) + y 3 2,C3l (x 1 2, " , «x 2 2fcal )] / Det Equation 6 

b = [(x t 2,ca, } 2 (y 2 2,ca, -y 3 2,c ") - (x 2 2,ca, ) 2 (y /^'.y^) + (x 3 2fca, ) 2 (y 1 2fca, -Y 2 2,C8, )] / Det Equation 7 

where 

Det - U^fix^x, 2 * - ) . {x 2 2fcal ) 2 (x 1 2,cal -x 3 2,ea, ) + (x 3 2,cal ) 2 (x 1 2,ca, .x 2 2,cal ) 

In block 216, the value of the predistortion parameter G q for the next iteration is determined as shown in Equation 8 
where the minimum of the polynomial y 2 ' cal occurs. 

GJn+1) - -b/(2a) Equation 8 

Note that the determinant (Det) cancels out of Equation 8. Therefore, the value of the determinant does not have to 
be calculated in an actual implementation. 

In block 218, the predistortion parameters b q and G q for the next iteration are used to collect the fourth, fifth 
and sixth set of samples. The fourth set of samples is collected with the value of the predistortion parameter K set 
equal to the previous value offset on the low side by the selected step value ~K. For notation purposes, we shall refer 
to this value as x 4 2,cal . The fourth set of samples is also collected with the predistortion parameter b { equal to the 
previous value offset on the low side by the selected step value Jb t . For notation purposes, we shall refer to this value 
as x 4 fcal . The fifth set of samples is collected by setting each of the predistortion parameters K and bj equal to their 
previous values. For notation purposes, we shall refer to the previous value of the predistortion parameter K(n) as 
x 5 2fMl and the previous value of the predistortion parameter bj(n) as x 5 ,tat . The sixth set of samples is collected with the 
predistortion parameter K set equal to its previous value offset on the high side by the selected step value ~K. The 
sixth set of samples is also collected by setting the predistortion parameter b ; equal to the previous value offset on the 
high side by the selected step value ~bj. For the ease of notation, we shall refer to these values as x 6 2iti] and x 6 ,cal , 
respectively. 

For the second three sets of samples, in block 220, the energy at f ca( is determined. For notation purposes, 
we shall refer to each of the energy values as y 4 fcal , y 5 fC8, f and y 6 fcat for sample sets 4, 5 and 6, respectively. In block 
222, three quadratic equations are developed as shown below in Equation 9. 
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y," - alx/" 1 ) 2 + bx/"' + c for i « 4, 5 and 6 Equation 9 

Also in block 222, the value of a and b are determined. For example, the formulas given below in Equations 10 and 1 1 
may be used to solve for the value of a and b based upon the three quadratic equations given by Equation 9. 

a = (y 4 !ca, (x 5 ,ca, -x 6 ,ca, ) - y 5 fca, (x 4 ,ca, -x 6 ,ca ') + ye'^x/"^^)] / Det Equation 10 



b- [(x 4 f " , ) z (y 5 fcaJ -y6 ,cal ) - (x 5 fca r (^V) + {x 6 ,ca, ) 2 (y 4 ' ca, .y 5 ,ca, )] / Det Equation 1 1 



where 



Det « (x 4 fca, ) 2 (x 5 fca, .x 6 ,ca, ) • (x 5 fca, ) 2 (x 4 ,ca, .x 6 , " , ) + (x 6 fcal ) 2 (x/ cal -x 5 fC31 ) 

In block 224, the value of the predistortion parameter b, for the next iteration is determined as shown in Equation 1 2 
where the minimum value of the polynomial y fcal occurs. 

10 Wn+1) «-b/(2a) Equation 12 

Note that the determinant (Det) cancels out of Equation 12. Therefore, the value of the determinant does not have to 
be calculated in an actual implementation. 

Using the second three sets of samples, in block 226, the energy at 2*f cal is determined. For notation 
purposes, we shall refer to the energy values as y 4 2fcal , y 5 2,cal , and y 6 2fcal for sample sets 4, 5 and 6, respectively. In 
15 block 228, three quadratic equations are developed as shown in Equation 13. 

y *ca. . a(x 2f«! )2 + bXj2 ™ + c for j = 4# 5 and 6 Equation 13 

Also, in block 228, the value of a and b are determined. For example, the formulas given in Equations 14 and 15 may 
be used to solve for the value of a and b based upon the three quadratic equations given by Equation 13. 

a = [y 4 2,C2, (x 5 2,cal -x 6 2fca, ) • y5 2fcal (x 4 2fca, .x 6 2,ca, ) + y 6 2fca, (x 4 2fcal -x 5 2,ca! )] / Det Equation 14 

20 b = [(x 4 2fca, ) 2 (y 5 2,cal -y 6 2,ca1 ) ■ (x 5 2 ' ca, ) 2 (y 4 2,cal .y 6 2,ca, ) + {x 6 2fca, ) 2 (y 4 2,cal -y 5 2f " ! )] / Det Equation 1 5 

where 

Det - (x 4 2fca, ) 2 (x 5 2,ca, -x s 2,cal ) - (x 5 2fca, ) 2 (x 4 2fcal -x 6 2fca ') + (x 6 2fcal ) 2 {x 4 2fca, •x 5 2fca, ) 

In block 228, the value of the crosstalk gain for the next iteration is determined as shown in Equation 16 where the 
minimum of the polynomial y 2fcal occurs. 

25 K(n+1) = -b/(2a) Equation 16 

Note that the determinant (Det) cancels out of Equation 16. Therefore, the value of the determinant does not have to 
be calculated in an actual implementation. 

In block 232, another set of samples is collected using the newly determined predistortion values, b q (n + 1), 
G q (n+1), qj(n + 1) and K(n+1). In block 232, the energy at f ca , and 2*f ca , corresponding to the verification samples is 
30 determined. Block 234 determines whether the ratio of the energy verification samples to the energy with no 
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predistortion exceeds a predetermined threshold. If so, the new predistortion values are deemed sufficient to be used 
and are applied during subsequent transmission mode operation. Flow terminates in block 238. If in block 234 it is 
found that the -ratio does not exceed the predetermined threshold, the flow process continues back to block 240 
where the value of n is incremented. In an alternative embodiment, the values of Ab q , AG q , Abj and AK are also 
5 changed, typically to smaller values, in block 240. From block 240, flow continues back to block 204 where a 

subsequent iteration is executed. 

In comparison with the prior art, the invention is robust because the quadratic polynomial minimization 
algorithm finds the value of the predistortion parameters which minimize the energy at the corresponding frequency. 
In contract, using a binary search method, if the initial estimate of the range of uncertainty does not include the value 
10 which corresponds to the minimum, the method does not return the exact value of the predistortion parameter which 

minimizes the energy. The invention typically converges within a single iteration while binary search methods may 
require a substantial number iterations.. 

Upon examination of the information herein, a myriad of alternative embodiments will be made readily 
available to one skilled in the art including the simple rearrangement of steps. For example, in one embodiment, the 
15 determination of the energy at f cal and 2 # f ca! may be accomplished simultaneously such that the processes shown in 

blocks 206 and 212 are performed at once. The parameters may be interchanged. For example, the predistortion 
parameter G can operate on the I channel. The grouping of the predistortion parameter can be varied such that the 
predistortion parameter b q is varied at the same time as the predistortion parameter K. In one embodiment, a 
waveform which comprises a sinusoidal component as well as other components may be applied during calibration 
20 mode. For example, a square-wave may be applied at frequency f„,. 

In one embodiment, the samples collected in block 232 of Figure 3 are the same some of the samples 
collected in block 204 during a subsequent iteration and the corresponding energy values are re-used rather than re- 
taken. The general principles of the invention may be combined with a variety of different convergence detection 
techniques. In Figure 3, convergence is detected by comparison of the spurious energies produced using the current 
25 values of the predistortion parameters with an initial spurious energy measurement. In alternative embodiments, 

convergence may be determined based upon the number of iterations which have been completed, the absolute level of 
the spurious energy, the change in the spurious energy between the last two iterations or a variety of other ways. 

Although the invention was described with respect to a direct conversion quadrature modulator, the 
teachings may be directly applied to other quadrature modulators such as those employing multiple conversion stages. 
30 In addition, the teachings of the invention may be extended to higher order spurious energy. 

The invention may be embodied in other specific forms without departing from its spirit or essential 
characteristics. The described embodiment is to be considered in all respects only as illustrative and not restrictive 
and the scope of the invention is, therefore, indicated by the appended claims rather than by the foregoing description. 
All changes which come within the meaning and range of equivalency of the claims are to be embraced within their 
35 scope. 
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WHAT IS CLAIMED IS : 

1 A method of determining a set of predistortion parameters for use with a quadrature modulator . 

comprising the steps of: 

(a) applying a sinusoidal wave to said quadrature modulator at a frequency f cal ; 
5 (b) collecting a first series of sets of samples reflecting changes to the value of a first set of predistortion 

parameters; 

(c) identifying a first series of energy values output by said quadrature modulator at said frequency f ca , and 
the frequency 2*f ca „ each of said series of energy values corresponding to one of said sets of samples; and 

(d) determining a subsequent value of said first set of predistortion parameters based upon a quadratic 
10 relationship between the value of said first set of predistortion parameters and said first series of energy values. 

2. The method of Claim 1 f further comprising the steps of: 

(e) collecting a second series of sets of samples reflecting changes to the value of a second set of 
predistortion parameters; 

15 (f) identifying a second series of energy values output by said quadrature modulator at said frequency f cal 

and said frequency 2*f cal , each of said series of energy values corresponding to one of said sets of samples; and 

(g) determining a subsequent value of said second set of predistortion parameters based upon a quadratic 
relationship between the value of said second set of predistortion parameters and said second series of energy values. 

20 3. The method of Claim 2, wherein the steps (b) ■ (g) are iteratively repeated until said first and 

second set of predistortion parameters converge. 

4. A method of determining a set of predistortion parameters for use in a predistortion block which 
drives an imperfect quadrature modulator in which a predistortion parameter, b ]f which determines a D.C. offset 

25 adjustment to an I channel, a predistortion parameter, b q , which determines a D.C. offset adjustment to an Q channel, 

a predistortion parameter, K, which determines a cross talk gain between said I channel and said Q channel, and a 
predistortion parameter, G, which determines a gain offset adjustment between said ! channel and said Q channel, the 
method comprising the steps of: 

applying a sinusoidal wave at frequency f ca , to said predistortion block; 

30 collecting a first set of samples of an output of said imperfect quadrature modulator corresponding to 

setting said predistortion parameter, G, to a lower value than a previous G value equal to x, 2 '" 1 , a first one of said 
predistortion parameters which determine said D.C. offset to a lower value than a previous b, value equal to x/ cal , 
setting said predistortion parameter, K, to a previous K value and a second one of said predistortion parameters which 
determines said D.C. offset to a previous b 2 value; 
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determining an energy level at the frequency f caI and the frequency 2f cal for said first set of samples equal to 
y/ cal and v, 21 * respectively; 

collecting a second set of samples corresponding to setting said predistortion parameter, G, to said previous 
G value equal to x 2 2,cal , said first one of said predistortion parameters which determine said D.C. offset to said previous 
b, value equal to x 2 ' cal , said predistortion parameter, K, to said previous K value and said second one of said 
predistortion parameters which determines said D.C. offset to said previous b 2 value; 

determining an energy level at the frequency f caI and the frequency 2f cal for said second set of samples equal 
to y 2 fcal and y 2 2fcal , respectively; 

collecting a third set of samples corresponding to setting said predistortion parameter, G, to a higher value 
than said previous G value equal to x 3 2fcal , said first one of said predistortion parameters which determine said D.C. 
offset to a higher value than said previous b, value equal to x 3 ,ca ', said predistortion parameter, K, to said previous K 
value and said second one of said predistortion parameters which determines said D.C. offset to said previous b 2 
value; 

determining an energy level at the frequency f cal and the frequency 2f cal for said third set of samples equal to 
y 3 fcaI and y 3 2,cal , respectively; 

determining a subsequent value of said predistortion parameter, G, based upon a quadratic relationship 

between said energy level y, 2 '" 1 and said value x, 2fcat , between said energy level y 2 2,caI and said value x 2 2,cal and 

between said energy level y 3 2fca! and said value x 3 2fcal ; 

determining a subsequent value of said first one of said predistortion parameters which determine said D.C. 

offset based upon a quadratic relationship between said energy level y/ caI and said value x/"' , between said energy 

level y 2 fcal and said value x 2 ' cal and between said energy level y 3 fcal and said value x 3 ,caI ; 

collecting a fourth set of samples of an output of said imperfect quadrature modulator corresponding to 

setting said predistortion parameter, K, to a lower value than a previous K value equal to x 4 2fC8 ', a second one of said 

predistortion parameters which determine said D.C. offset to a lower value than a previous b 2 value equal to x 4 caI , 

setting said predistortion parameter, G, to said subsequent value thereof and said first one of said predistortion 

parameters which determines said D.C. offset to said subsequent value thereof- 
determining an energy level at the frequency f caJ and the frequency 2f cal for said fourth set of samples equal 

to y 4 fcal and y 4 2fcaf , respectively; 

collecting a fifth set of samples corresponding to setting predistortion parameter, K, to said previous K value 

equal to x 4 2,cal , said second one of said predistortion parameters which determine said D.C. offset to said previous b 2 

value equal to x 4 "', said predistortion parameter, G, to said subsequent value thereof and said first one of said 

predistortion parameters which determines said D.C. offset to said subsequent value thereof- 
determining an energy level at the frequency f caJ and the frequency 2f caI for said fifth set of samples equal to 

y 5 caI and y 5 2fcal , respectively; 
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collecting a sixth set of samples corresponding to setting said predistortion parameter, K, to a higher value 
than said previous K value equal to x 6 2,cal , said second one of said predistortion parameters which determine said D.C. - 
offset to a higher value than said previous b 2 value equal to x 6 fcal , said predistortion parameter, G, to said subsequent 
value thereof and said first one of said predistortion parameters which determines said D.C. offset to said subsequent 
5 value thereof; 

determining an energy level at the frequency f ca , and the frequency 2f cal for said sixth set of samples equal to 
y 6 fcaI and y 6 2fcal , respectively; 

determining a subsequent value of said predistortion parameter, K, based upon a quadratic relationship 
between said energy level y 4 2fcaI and said value x 4 2fcal , between said energy level y 5 2fcaI and said value x 5 2fcal and 
10 between said energy level y 6 2,cal and said value x 6 2fca1 ; and 

determining a subsequent value of said second one of said predistortion parameters which determine said 
D.C. offset based upon a quadratic relationship between said energy level y 4 ca> and said value x 4 cal , between said 
energy level y 5 cal and said value x 5 cal and between said energy level y 6 fcal and said value x 6 fcal . 

15 5. A method of determining a set of predistortion parameters for use in predistortion block which 

drives an imperfect quadrature modulator in which a predistortion parameter, b i# which determines a D.C. offset 
adjustment to an I channel, a predistortion parameter, b q , which determines a D.C. offset adjustment to a Q channel, a 
predistortion parameter, K, which determines a cross talk gain between said I channel and said Q channel, and a 
predistortion parameter, G, which determines a gain offset adjustment between said I channel and said Q channel, said 

20 method comprising the steps of: 

applying a sinusoidal wave at frequency f ca , to said predistortion block; 

setting said predistortion parameter, G, equal to a previous value, G(n), offset on the low side by the 
selected step value ~G wherein G(n) • AG - x, 2fcal ; 

setting said predistortion parameter, b q , equal to a previous value, b q (n), offset on the low side by the 
25 selected step value ~b q wherein b q (n) - Ab q - x/ cal ; 

setting said predistortion parameter, K, equal to a previous value, K(n); 
setting said predistortion parameter, b jr equal to a previous value, b ; (n); 

collecting a first set of samples of an output of said imperfect quadrature modulator and determining an 
energy level at the frequency f ca , and the frequency 2f cal for said first set of samples equal to y/ caI and y, 2fcal , 
30 respectively; 

setting said predistortion parameter, G, equal to said previous value, G(n) - x 2 2,cal ; 
setting said predistortion parameter, b q , equal to said previous value, b q (n) - x 2 fcal ; 
collecting a second set of samples of said output of said imperfect quadrature modulator and determining an 
energy level at the frequency f ca , and the frequency 2f caI for said second set of samples equal to y 2 ,cal and y 2 2fcal , 
35 respectively; 
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setting said predistortion parameter, G, equal to said previous value, G(n), offset on the high side by the 
selected step value ~G wherein G(n) + AG - x 3 2fcal ; 

setting said predistortion parameter, b q , equal to said previous value, b q {n), offset on the high side by the 
selected step value ~b q wherein b q (n) + Ab q - x 3 ,caJ ; 

collecting a third set of samples of said output of said imperfect quadrature modulator and determining an 
energy level at the frequency f cal and the frequency 2f cal for said third set of samples equal to y 3 ,cal and y 3 2,M, # 
respectively; 

determining a subsequent value of said predistortion parameter, b q , equal to b q (n + 1) by solving 

b q (n+1)--b/{2a) 

where 

y tea. _ a{x fc3i )2 + bx fca. + c for J . 1# 2 and 3; 

determining a subsequent value of said predistortion parameter, G, equal to G(n + 1) by solving 

G(n + 1)--b/{2a) 

where 

y IW m a(x 2fca. )2 + bx 2fcal + c f Qr j _ 2 and 3 

setting said predistortion parameter, G, equal to said next value, G(n + In- 
setting said predistortion parameter, b q , equal to said next value, b q (n + In- 
setting said predistortion parameter, K, equal to said previous value, K(n) offset on the low side by the 

selected step value _K wherein K(n) - AK - x 4 2feal ; 

setting said predistortion parameter, b u equal to said previous value, bj(n) offset on the low side by the 

selected step value J>j wherein b;(n) • Abj = x 4 fcal ; 

collecting a fourth set of samples of said output of said imperfect quadrature modulator and determining an 

energy level at the frequency and the frequency 2f cal for said fourth set of samples equal to y 4 fCfll and y 4 2fcal , 

respectively; 

setting said predistortion parameter, K, equal to said previous value, K(n) - x 5 2fcal ; 
setting said predistortion parameter, bj, equal to said previous value, b,(n) - x 5 ' cal ; 

collecting a fifth set of samples of said output of said imperfect quadrature modulator and determining an 
energy level at the frequency f ca , and the frequency 2f taJ for said fifth set of samples equal to y 5 ,cal and y 5 2fcaI , 
respectively; 

setting said predistortion parameter, K, equal to said previous value, K(n), offset on the high side by the 
selected step value ~K wherein K(n) + AK - x 6 2fcal ; 
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setting said predistortion parameter, b t , equal to said previous value, b;(n), offset on the high side by the 
selected step value ~bj wherein b^n) + Ab s - x 6 fcal ; 

""collecting a sixth set of samples of said output of said imperfect quadrature modulator and determining an 
energy level at the frequency f caJ and the frequency 2f caI for said sixth set of samples equal to y 6 ,cal and y 6 2fcaI , 
respectively; 

determining a subsequent value of said predistortion parameter, b ir equal to bj(n + 1) by solving 

Mn+D-.b/tfa) 

where 

y « gfx/" 1 ) 2 + bx^ 1 + c for i « 4, 5 and 6; and 
determining a subsequent value of said predistortion parameter, K, equal to K(n + 1) by solving 

K(n+1) = -b/(2a) 

where 

y 2fca« _ g(x 2fcal,2 + ^fca! + c f or j = 4 , 5 gnd B- 

6. The method of determining said set of predistortion parameters of Claim 5, wherein said imperfect 
quadrature modulator is a direct conversion modulator. 

7. The method of determining said set of predistortion parameters of Claim 5 r wherein said step of 
applying a sinusoidal wave at frequency f ca , to said predistortion block and said steps of collecting said first, second, 
third, fourth, fifth and six sample sets occur during a calibration mode and wherein said steps of determining said 
energy levels at the frequency f M( and the frequency 2f cal and said determining said subsequent values occur during 
said calibration mode or during a transmit mode. 

8. The method of determining said set of predistortion parameters of Claim 5, further comprising the 

steps of: 

setting said predistortion parameter, G, equal to said next value, G(n+ In- 
setting said predistortion parameter, b q , equal to said next value, b q (n + In- 
setting said predistortion parameter, K, equal to said next value, K{n+ In- 
setting said predistortion parameter, b [f equal to said next value, b,(n + 1); 

collecting a seventh set of samples of said output of said imperfect quadrature modulator and determining an 
energy level at the frequency f ca , and the frequency 2f cat ; 

comparing said energy values at the frequency f ca1 and the frequency 2f cal to previously collected energy 
values at the frequency f ca , and the frequency 2f cal ; and 
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setting said previous values equal to said subsequent values and repeating said steps of collecting samples 
and determining said subsequent values if a difference between said energy values and said previously collected set of 
energy values is less than a predetermined value. 

9. An apparatus for determining a set of predistortion parameters for use in a predistortion block 
which drives an imperfect quadrature modulator in which a predistortion parameter, b lt which determines a D.C. offset 
adjustment to an I channel, a predistortion parameter, b q , which determines a D.C. offset adjustment to an Q channel, 
a predistortion parameter, K, which determines a cross talk gain between said I channel and said Q channel, and a 
predistortion parameter, G, which determines a gain offset adjustment between said I channel and said Q channel, the 
apparatus comprising: 

means for applying a sinusoidal wave at frequency f ca , to said predistortion block; 

means for collecting a first set of samples of an output of said imperfect quadrature modulator 
corresponding to setting said predistortion parameter, G, to a lower value than a previous G value equal to x, 2 '" 1 , a 
first one of said predistortion parameters which determine said D.C. offset to a lower value than a previous b, value 
equal to x/ ca, f setting said predistortion parameter, K, to a previous K value and a second one of said predistortion 
parameters which determines said D.C. offset to a previous b 2 value; 

means for determining an energy level at the frequency f cal and the frequency 2f ca , for said first set of 
samples equal to y/ cal and y, 2 '" 1 , respectively; 

means for collecting a second set of samples corresponding to setting said predistortion parameter, G, to 
said previous G value equal to x 2 2fcal , said first one of said predistortion parameters which determine said D.C. offset to 
said previous b 1 value equal to x 2 fcaI , said predistortion parameter, K, to said previous K value and said second one of 
said predistortion parameters which determines said D.C. offset to said previous b 2 value- 
means for determining an energy level at the frequency f ca , and the frequency 2f cal for said second set of 
samples equal to y 2 ,cal and y 2 2fcal , respectively; 

means for collecting a third set of samples corresponding to setting said predistortion parameter, G, to a 
higher value than said previous G value equal to x 3 2fcal , said first one of said predistortion parameters which determine 
said D.C. offset to a higher value than said previous b, value equal to x 3 fcal , said predistortion parameter, K, to said 
previous K value and said second one of said predistortion parameters which determines said D.C. offset to said 
previous b 2 value; 

means for determining an energy level at the frequency f C3f and the frequency 2^ for said third set of 
samples equal to y 3 ,cat and y 3 2fcal , respectively; 

means for determining a subsequent value of said predistortion parameter, G, based upon a quadratic 
relationship between said energy level y, 2,ca} and said value x, 21 " 1 , between said energy level y 2 2,cal and said value x 2 2,cal 
and between said energy level y 3 2,cal and said value x 3 2,cal ; 
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means for determining a subsequent value of said first one of said predistortion parameters which determine 
said D.C. offset based upon a quadratic relationship between said energy level y/ ca1 and said value x/ ca1 , between said - 
energy level y 2 fcal and said value x 2 fcal and between said energy level y 3 fcal and said value x 3 fcal ; 

means for collecting a fourth set of samples of an output of said imperfect quadrature modulator 
5 corresponding to setting said predistortion parameter, K, to a lower value than a previous K value equal to x 4 2fcal , a 

second one of said predistortion parameters which determine said D.C. offset to a lower value than a previous b 2 value 
equal to x 4 cai , setting said predistortion parameter, G, to said subsequent value thereof and said first one of said 
predistortion parameters which determines said O.C. offset to said subsequent value thereof; 

means for determining an energy level at the frequency f ca , and the frequency 2f ca , for said fourth set of 
1 0 samples equal to y 4 fca) and y 4 2fcal , respectively; 

means for collecting a fifth set of samples corresponding to setting predistortion parameter, K, to said 
previous K value equal to x 4 2fcal , said second one of said predistortion parameters which determine said O.C. offset to 
said previous b 2 value equal to x 4 ca! , said predistortion parameter, G, to said subsequent value thereof and said first one 
of said predistortion parameters which determines said D.C. offset to said subsequent value thereof; 
1 5 means for determining an energy level at the frequency f cal and the frequency 2f ca , for said fifth set of 

samples equal to y 5 cal and y 5 2fcaJ , respectively; 

means for collecting a sixth set of samples corresponding to setting said predistortion parameter, K, to a 
higher value than said previous K value equal to x 6 2fca> , said second one of said predistortion parameters which 
determine said D.C. offset to a higher value than said previous b 2 value equal to x 6 fcaJ , said predistortion parameter, G, 
20 to said subsequent value thereof and said first one of said predistortion parameters which determines said D.C. offset 

to said subsequent value thereof; 

determining an energy level at the frequency f cal and the frequency 2f ca , for said sixth set of samples equal to 
y 8 fcal and y 6 2fcal , respectively; 

determining a subsequent value of said predistortion parameter, K, based upon a quadratic relationship 
25 between said energy level y 4 2,cat and said value x 4 2fcaJ , between said energy level y 5 2fcal and said value x 5 2fcal and 

between said energy level y 6 2fcal and said value x 6 2fcal ; and 

determining a subsequent value of said second one of said predistortion parameters which determine said 
D.C. offset based upon a quadratic relationship between said energy level y 4 cal and said value x 4 cal , between said 
energy level y 5 cal and said value x 5 cal and between said energy level y 6 ,tal and said value x 6 fcal . 

30 

10. A transmit modulator comprising: 

a predistortion block producing a predistorted output signal comprising: 

a D.C. offset compensation summer which determines a D.C. offset adjustment to an I 
channel based upon a predistortion parameter, b t ; 
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a D.C. offset compensation summer which determines a D.C. offset adjustment to an Q 
channel based upon a predistortion parameter, b q ; 

a cross talk amplifier which determines a cross talk gain between said I channel and said Q 
channel based upon a predistortion parameter, K; and 

a gain compensation amplifier which determines a gain offset adjustment between said I 
channel and said Q channel based upon a predistortion parameter, G; 

an imperfect quadrature modulator coupled to said predistortion block and receiving said predistorted output 
signal and producing a modulated signal; 

a transformer coupled to said imperfect quadrature modulator and receiving said modulated signal and 
producing a digital representation of a spectrum of said modulated signal; 

a quadratic polynomial minimization calculator which determines a subsequent value of said predistortion 
parameters, b u b q , K and G based upon a quadratic relationship between energies present in said digital representation 
of said spectrum of said modulated signal and the values of said predistortion parameter b ir b q , K and G. 

11. The transmit modulator of Claim 10 wherein said quadratic polynomial minimization calculator 
contains a series of processes that when executed perform the steps of: 

collecting a first set of samples of an output of said imperfect quadrature modulator corresponding to 
setting said predistortion parameter, G, to a lower value than a previous G value equal to x, 2 '" 1 , a first one of said 
predistortion parameters which determine a D.C. offset to a lower value than a previous b, value equal to x/ Ml , setting 
said predistortion parameter, K, to a previous K value and a second one of said predistortion parameters which 
determines a D.C. offset to a previous b 2 value; 

determining an energy level at the frequency f caf and the frequency 2f cal for said first set of samples equal to 
y/ ca ' and y, 2 '" 1 , respectively; 

collecting a second set of samples corresponding to setting said predistortion parameter, G, to said previous 
G value equal to x 2 2fC3, f said first one of said predistortion parameters which determine said D.C. offset to said previous 
b, value equal to x 2 ,cal , said predistortion parameter, K, to said previous K value and said second one of said 
predistortion parameters which determines said D.C. offset to said previous b 2 value; 

determining an energy level at the frequency f caJ and the frequency 2f cal for said second set of samples equal 
to y 2 ,cal and y 2 2,cal , respectively; 

collecting a third set of samples corresponding to setting said predistortion parameter, G, to a higher value 
than said previous G value equal to x 3 2fcal , said first one of said predistortion parameters which determine said D.C. 
offset to a higher value than said previous b, value equal to x 3 fca! , said predistortion parameter, K, to said previous K 
value and said second one of said predistortion parameters which determines a D.C. offset to said previous b 2 value; 

determining an energy level at the frequency f caI and the frequency 2f cal for said third set of samples equal to 
y 3 fcaf and y 3 2,ca ', respectively; 
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determining a subsequent value of said predistortion parameter, G, based upon a quadratic relationship 

between said energy level y, 2fcal and said value x* Ci] , between said energy level y 2 2lcal and said value x 2 2fcal and 

between said energy level y 3 2fcal and said value x 3 2fcal ; 

determining a subsequent value of said first one of said predistortion parameters which determine said D.C. 
5 offset based upon a quadratic relationship between said energy level y/ EaI and said value x/ cal , between said energy 

level y 2 ,cal and said value x 2 ,cal and between said energy level y 3 ,cal and said value x 3 ,cal ; 

collecting a fourth set of samples of an output of said imperfect quadrature modulator corresponding to 

setting said predistortion parameter, K, to a lower value than a previous K value equal to x 4 2fcal , a second one of said 

predistortion parameters which determine said D.C. offset to a lower value than a previous b 2 value equal to x 4 cal , 
10 setting said predistortion parameter, G, to said subsequent value thereof and said first one of said predistortion 

parameters which determines a D.C. offset to said subsequent value thereof- 
determining an energy level at the frequency f cal and the frequency 2f ca , for said fourth set of samples equal 

to y 4 ,caJ and y 4 2fca, r respectively; 

collecting a fifth set of samples corresponding to setting predistortion parameter, K, to said previous K value 
1 5 equal to x 4 2,cal , said second one of said predistortion parameters which determine said D.C. offset to said previous b 2 

value equal to x 4 ca, r said predistortion parameter, G, to said subsequent value thereof and said first one of said 

predistortion parameters which determines a D.C. offset to said subsequent value thereof; 

determining an energy level at the frequency f ca , and the frequency 2f M , for said fifth set of samples equal to 

y 5 " ! and y 5 2fca! f respectively; 

20 collecting a sixth set of samples corresponding to setting said predistortion parameter, K, to a higher value 

than said previous K value equal to x 6 2,caI , said second one of said predistortion parameters which determine said D.C. 
offset to a higher value than said previous b 2 value equal to x 6 ,cal , said predistortion parameter, G, to said subsequent 
value thereof and said first one of said predistortion parameters which determines a D.C. offset to said subsequent 
value thereof; 

25 determining an energy level at the frequency f cal and the frequency 2f ca , for said sixth set of samples equal to 

Ye™ and y 6 2,cal , respectively; 

determining a subsequent value of said predistortion parameter, K, based upon a quadratic relationship 

between said energy level y 4 2fcal and said value x 4 2fcal , between said energy level y 5 2!cal and said value x 5 2,cal and 

between said energy level y 6 2fcat and said value x 6 2,ca1 ; and 
30 determining a subsequent value of said second one of said predistortion parameters which determine said 

D.C. offset based upon a quadratic relationship between said energy level y 4 cal and said value x 4 ca> , between said 

energy level y 5 cal and said value x 5 cal and between said energy level y 6 ,cal and said value x 6 fcal . 

12. A programmed storage device storing a series of process that when executed perform a method of 
35 determining a set of predistortion parameters for use in a predistortion block which drives an imperfect quadrature 
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modulator in which a predistortion parameter, b j( which determines a D.C. offset adjustment to an I channel, a 
predistortion parameter, b,, which determines a D.C. offset adjustment to an Q channel, a predistortion parameter, K, 
which determines a cross talk gain between said I channel and said Q channel, and a predistortion parameter. G, which 
determines a gain offset adjustment between said I channel and said Q channel, the method comprising the steps of: 
applying a sinusoidal wave at frequency f,,, to said predistortion block; 

collecting a first set of samples of an output of said imperfect quadrature modulator corresponding to 
setting said predistortion parameter, G, to a lower value than a previous G value equal to x, 2 '" 1 , a first one of said 
predistortion parameters which determine said D.C. offset to a lower value than a previous b, value equal to x,* - , 
setting said predistortion parameter. K, to a previous K value and a second one of said predistortion parameters which 
determines said D.C. offset to a previous b 2 value; 

determining an energy level at the frequency f cal and the frequency 2f c „ for said first set of samples equal to 
y,"* 1 and y, 2 '"', respectively; 

collecting a second set of samples corresponding to setting said predistortion parameter, G, to said previous 
G value equal to x 2 2,Ml , said first one of said predistortion parameters which determine said D.C. offset to said previous 
b, value equal to x 2 M , said predistortion parameter, K, to said previous K value and said second one of said 
predistortion parameters which determines said D.C. offset to said previous b 2 value; 

determining an energy level at the frequency f Ml and the frequency 2f Ml for said second set of samples equal 
to y 2 '"' and y 2 2,ul , respectively; 

collecting a third set of samples corresponding to setting said predistortion parameter, G, to a higher value 
than said previous G value equal to x 3 2 "", said first one of said predistortion parameters which determine said D.C. 
offset to a higher value than said previous b, value equal to x 3 ,Ml , said predistortion parameter, K, to said previous K 
value and said second one of said predistortion parameters which determines said D.C. offset to said previous b 2 
value; 

determining an energy level at the frequency f eil and the frequency 2f c * for said third set of samples equal to 
y 3 "" and y^""', respectively; 

determining a subsequent value of said predistortion parameter, G, based upon a quadratic relationship 
between said energy level y, 2 "" and said value x, 2 '" 1 , between said energy level y, 2 """ and said value x, 2 *" and 
between said energy level y 3 2,c " and said value x 3 2 ""; 

determining a subsequent value of said first one of said predistortion parameters which determine said D.C. 
offset based upon a quadratic relationship between said energy level y," and said value x, ,M * , between said energy 
level y 2 lcal and said value x 2 ' cal and between said energy level y 3 '"' and said value x 3 ,tal ; 

collecting a fourth set of samples of an output of said imperfect quadrature modulator corresponding to 
setting said predistortion parameter, K, to a lower value than a previous K value equal to x 4 2 ' Ml , a second one of said 
predistortion parameters which determine said D.C. offset to a lower value than a previous b 2 value equal to x 4 eal , 
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setting said predistortion parameter, G, to said subsequent value thereof and said first one of said predistortion 
parameters which determines said D.C. offset to said subsequent value thereof; 

determining an energy level at the frequency f ca , and the frequency 2f cal for said fourth set of samples equal 
to y 4 fca! and y 4 2fcal , respectively; 

collecting a fifth set of samples corresponding to setting predistortion parameter, K, to said previous K value 
equal to x 4 2fcal , said second one of said predistortion parameters which determine said D.C. offset to said previous b 2 
value equal to x 4 ca, r said predistortion parameter, G, to said subsequent value thereof and said first one of said 
predistortion parameters which determines said D.C. offset to said subsequent value thereof; 

determining an energy level at the frequency f cal and the frequency 2f ca , for said fifth set of samples equal to 
y 5 caI and y 5 2fcal , respectively; 

collecting a sixth set of samples corresponding to setting said predistortion parameter, K, to a higher value 
than said previous K value equal to x 6 2fcal , said second one of said predistortion parameters which determine said D.C. 
offset to a higher value than said previous b 2 value equal to x 6 ,cal , said predistortion parameter, G, to said subsequent 
value thereof and said first one of said predistortion parameters which determines said D.C. offset to said subsequent 
value thereof; 

determining an energy level at the frequency f ca , and the frequency 2f cal for said sixth set of samples equal to 
y 6 ,taI and y 6 2kal , respectively; 

determining a subsequent value of said predistortion parameter, K, based upon a quadratic relationship 
between said energy level y 4 2fcaI and said value x 4 2,C8j , between said energy level y 5 2fcal and said value x 5 2fcal and 
between said energy level y 6 2fca ' and said value x 6 2fcal ; and 

determining a subsequent value of said second one of said predistortion parameters which determine said 
D.C. offset based upon a quadratic relationship between said energy level y 4 " ! and said value x 4 cal , between said 
energy level y 5 ca! and said value x s caI and between said energy level y 6 fcal and said value x 6 fcal . 



•24- 



WO 00/31941 



PCT/US99/27256 



1/3 



-to 



-/4 



-/6 



TM 




QM 










PA 



-/2 



*/4 



-/6 



W 



TM 




PD 




QM 




PA 













SUBSTITUTE SHEET (RULE 26) 



WO 00/31941 



2/3 



PCT/US99/27256 




SUBSTITUTE SHEET (RULE 26) 



WO 00/31941 



,-202 



START 



3/3 



^200 



DETERMINE THE ENERGY AT 
fcal AND 2f ca | WITH 
NO PREDISTORTION 



-204 



1 



COLLECT SETS OF SAMPLES AS 
SHOWN IN TABLE I 



-20S 



I 



DETERMINE THE ENERGY AT f ca | 
FOR SAMPLE SETS 1,2,3 = 
yi fcal .y 2 fcal AND y 3 fcal 



S-20S 



DETERMINE a AND 


b WHERE y; fc a' 


= a(x{fcal) 2 + 


bx;fcal + c 


for I = 


1.2,3 


S-2/0 i 



DETERMINE 
b q (n+1) = -b/(2a) 



^-2/2 



I 



DETERMINE THE ENERGY AT 2*f ca! 
FOR SAMPLE SETS 1,2,3 = y 1 2f cal 
y 2 2fcal AND y 3 2fcal 



DETERMINE a AND b WHERE 
^ 2fcal = a(x. 2fc al)2 + bxj 2fcal 
+ c, for i = 1,2,3 



-2/6 j 



DETERMINE 
G q (n+1) = -b/(2a) 



^-2/8 



I 



COLLECT SET OF SAMPLES AS 
SHOWN IN TABLE II 

DETERMINE THE ENERGY AT f ca | 

FOR SAMPLE SETS 4,5,6 = 
y 4 fcQl,y 5 feol AND ygfcal 



PCT/US99/27256 



1 



DETERMINE a AND b WHERE y.fcal 
= a(x.fcal)2 + bx .fcal + ' c 

for i = 4,5,6 



DETERMINE 
b,(n+1) = -b/(2a) 



DETERMINE THE ENERGY AT 2*f ca , 
FOR SAMPLE SETS 4,5,6 
y 4 2fcal,y 5 2fcal AND yg2fcal 



I 



-228 



DETERMINE a AND b WHERE 
y.2fcal= a (xj2fcal)2 + bx.2fcal 

+ c for i = 4,5,6 

I S-2JO 240 



DETERMINE 
K(n+1) = -b/(2a) 



j ^2J2 



n=n+1 



COLLECT A SET OF VERIFICATION 

SAMPLES USING 
Gq(n+1), b q (n+1) AND bj(n+1) 



DETERMINE THE ENERGY AT f ca . 
AND 2f ca | FOR THE 
VERIFICATION SAMPLES 




DOES THE RATIO OF THE 
ENERGY IN THE VERIFICATION 
SAMPLE TO THE ENERGY 
WITH NO PREDISTORTION 
EXCEED A PREDETERMINED 
THRESHOLD? 



j s- 238 




END 



r/aj 



SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



Int tional Application No 

PCT/US 99/27256 



A CLASSIFICATION OF SUBJECT MATTER , 

IPC 7 H04L27/36 H03C3/40 



According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 H04L H03C 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ° Citation of document, with indication, where appropriate, of the relevant passages 



LOHTIA A ET AL: "ADAPTIVE DIGITAL 
LINEARIZATION OF RF POWER" 
CANADIAN JOURNAL OF ELECTRICAL AND 
COMPUTER ENGINEERING, CANADIAN SOCIETY FOR 
ELECTRICAL ENGINEERING, 

vol. 20, no. 2, 1 April 1995 (1995-04-01), 
pages 65-71, XP000671568 
Victoria, BC, CANADA 

ISSN: 0840-8688 
page 67, left-hand column, line 16 
-right-hand column, line 18 

figure 4 



-/-- 



Relevant to claim No. 



1-3 



4-6,9, 
10,12 



| X [ frurtn9r documents are listed in the continuation of box C. 



Patent family members are listed in annex. 



° Special categories of cited documents : 

"A" document defining the general state of the art which is not 
considered to be of particular relevance 

"E H earlier document but published on or after the international 
filing date 

"L" document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



"T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

"&° document member of the same patent family 



Date of the actual completion of the international search 

22 March 2000 


Date of mailing of the international search report 

29/03/2000 


Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL-2280 HVRijswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 


Authorized officer 

Farese, L 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



tm itional Application No 

PCT/US 99/27256 



C.(Contlnuation) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category J 



Citation of document, with indication. where appropriate, of the relevant passages 



RolQvant to claim No. 



A 
A 



STAPLETON S P ET AL: "AN ADAPTIVE 

PREDISTORTER FOR A POWER AMPLIFIER BASED 

ON ADJACENT CHANNEL EMISSIONS" 

IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY 

IEEE INC. , 

vol. 41, no. 1, 

1 February 1992 (1992-02-01), pages 49-56, 

XP000297040 

New York, US 

ISSN: 0018-9545 
page 52, left-hand column, line 1 - line 
25 

page 53, right-hand column, line 1 - line 
5 

page 54, left-hand column, line 22 - line 
28 

figures 1,2 

HILBORN D S ET AL: "AN ADAPTIVE DIRECT 
CONVERSION TRANSMITTER" 

IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY 
IEEE INC. , 

vol. 43, no. 2, 1 May 1994 (1994-05-01), 
pages 223-233, XP000460100 
New York, US 

ISSN: 0018-9545 
page 228, right-hand column, line 39 - 
line 46 

page 229, left-hand column, line 31 - line 
49 

US 5 524 285 A (WRAY ANTHONY J ET AL) 
4 June 1996 (1996-06-04) 
column 5, line 42 - line 53 
figure 2B 



1-3 



1-3 



4-6,9, 
10,12 

1 



Four PCT/ISA/21 0 (continuation ol second sheet) (July 1 892) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

Information on patent family members 



Int jonal Application No 

PCT/US 99/27256 



Patent document 


Publication 


Patent family 


Publication 


cited in search report 


date 


member(s) 


date 


US 5524285 ... A 


04-06-1996 


NONE 





Foim PCT/ISA/210 (patent family annex) (July 1992) 



